*  ★ 


NFS 

PRAlSTANTlAPERSCIENm^ 


NAVAL 

POSTGRADUATE 

SCHOOL 


MONTEREY,  CALIFORNIA 

THESIS 


EFFECTS  OF  A  SUSPENDED  SEDIMENT  LAYER  ON  ACOUSTIC 

IMAGERY 

by 

Michael  Cornelius 

June  2004 

Thesis  Advisor: 
Second  Reader: 

Peter  C.  Chu 

Melvin  Wagstaff  (NAVO) 

Approved  for  public  release;  distribution  is  unlimited 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


REPORT  DOCUMENTATION  PAGE 


Form  Approved  0MB  No. 


0704-0188 


Public  reporting  burden  for  this  colleetion  of  information  is  ( 
response,  ineluding  the  time  for  reviewing  instruetion,  searehin^ 
and  maintaining  the  data  needed,  and  eompleting  and  reviewii 
Send  eomments  regarding  this  burden  estimate  or  any  oth 
information,  ineluding  suggestions  for  reducing  this  burden,  to  W 
Direetorate  for  Information  Operations  and  Reports,  1215  Jeffei 
Arlington,  VA  22202-4302,  and  to  the  Offiee  of  Management  ai 
Project  (0704-0188)  Washington  DC  20503. 

estimated  to  average  1  hour  per 
;  existing  data  sources,  gathering 
ig  the  collection  of  information, 
er  aspect  of  this  collection  of 
Washington  headquarters  Services, 
•son  Davis  Highway,  Suite  1204, 
ad  Budget,  Paperwork  Reduction 

1.  AGENCY  USE  ONLY  (Leave  blank) 

2.  REPORT  DATE 

June  2004 

3.  REPORT  TYPE  AND  DATES 
COVERED 

Master' s  Thesis 

4.  TITLE  AND  SUBTITLE:  Effects  of  a  Suspended  Sediment  Layer  on 
Acoustic  Imagery 

5.  FUNDING  NUMBERS 

6.  AUTHOR(S)  Cornelius,  Michael 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Postgraduate  School 

Monterey,  CA  93943-5000 

8.  PERFORMING 

ORGANIZATION  REPORT 
NUMBER 

9.  SPONSORING  /MONITORING  AGENCY  NAME(S)  AND 
ADDRESS(ES) 

NAVOCEANO 

Stennis  Space  Center,  MS  39522-5001 

10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 

11.  SUPPLEMENTARY  NOTES  The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect  the 
official  policy  or  position  of  the  Department  of  Defense  or  the  U.S.  Government. 

12a.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited 

12b.  DISTRIBUTION  CODE 

13.  ABSTRACT  (maximum  200  words) 

The  Navy’s  CASS/GRAB  sonar  model  is  used  to  accurately  simulate  a  side-scan  sonar  image  with 
a  mine-like  object  present  through  its  reverberation  characteristics.  The  acoustic  impact  of  a  suspended 
sediment  layer  is  investigated  numerically  using  CASS/GRAB  through  changing  the  volume  scattering 
characteristics  of  the  lower  water  column.  A  range  of  critical  values  of  volume  scattering  strength  were 
discovered  through  repeated  model  simulations.  An  understanding  of  the  acoustic  characteristics  of 
suspended  sediment  layers  can  aid  the  Navy  in  the  detection  of  mines  that  might  exist  within  these  layers. 

14.  SUBJECT  TERMS 

CASS/GRAB,  Modeling  and  Simulation,  Oceanography,  Mine  Warfare 

Nepheloid  Layer,  Bottom  Boundary  Layer 

15.  NUMBER  OF 
PAGES  64 

16.  PRICE  CODE 

17.  SECURITY 
CLASSIFICATION  OF 
REPORT 

Unclassified 

18.  SECURITY 
CLASSIFICATION  OF  THIS 
PAGE 

Unclassified 

19.  SECURITY 
CLASSIFICATION 

OF  ABSTRACT 

Unclassified 

20.  LIMITATION 

OF  ABSTRACT 

UL 

NSN  7540-01-280-5500  Standard  Form  298  (Rev.  2-89) 


Prescribed  by  ANSI  Std.  239- 1 8 


1 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


11 


Approved  for  public  release;  distribution  is  unlimited 


EFFECTS  OF  A  SUSPENDED  SEDIMENT  LAYER  ON  ACOUSTIC  IMAGERY 

Michael  Cornelius 
Lieutenant,  United  States  Navy 
B.S.,  United  States  Naval  Academy,  1996 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  METEOROLOGY  AND  PHYSICAL 

OCEANOGRAPHY 


from  the 


NAVAL  POSTGRADUATE  SCHOOL 
June  2004 


Author:  Michael  Cornelius 


Approved  by:  Peter  C.  Chu 

Thesis  Advisor 


Melvin  D.  Wagstaff  (NAVO) 
Second  Reader 


Mary  Batteen 

Chairman,  Department  of  Oceanography 


111 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


IV 


ABSTRACT 


The  Navy’s  CASS/GRAB  sonar  model  is  used  to  accurately  simulate  a  side-scan 
sonar  image  with  a  mine-like  object  present  through  its  reverberation  characteristics.  The 
acoustic  impact  of  a  suspended  sediment  layer  is  investigated  numerically  using 
CASS/GRAB  through  changing  the  volume  scattering  characteristics  of  the  lower  water 
column.  A  range  of  critical  values  of  volume  scattering  strength  were  discovered  through 
repeated  model  simulations.  An  understanding  of  the  acoustic  characteristics  of 
suspended  sediment  layers  can  aid  the  Navy  in  the  detection  of  mines  that  might  exist 
within  these  layers. 
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I. 


INTRODUCTION 


In  the  fifteen  years  since  the  fall  of  the  Berlin  Wall  and  the  end  of  the  “Cold 
War”,  United  States  Naval  Forces  are  still  transitioning  from  a  large  “blue  water”  navy  to 
a  smaller,  more  mobile  force  capable  of  combat  in  the  littoral  regions  of  the  world.  The 
looming  specter  of  the  Soviet  Union  has  been  replaced  by  smaller  regional  threats  that 
utilize  unconventional  tactics  to  inflict  damage  on  those  they  oppose. 

One  such  tactic  that  was  prevalent  during  Operation  Iraqi  Freedom  was  that  of 
using  sea  mines  to  block  waterways.  Naval  mines  are  inexpensive,  easily  deployed,  and 
possess  extraordinary  destructive  capability.  Once  deployed,  mine  fields  require  very 
little  maintenance  and  can  effectively  block  access  to  ports,  deter  amphibious  operations, 
prevent  humanitarian  assistance  and  psychologically  affect  those  assigned  to  patrol  the 
waters  in  the  vicinity  of  the  minefield.  The  damage  caused  by  naval  mines  to  U.S.  assets 
is  well  documented.  The  USS  SAMUEL  B.  ROBERTS  (FEG-59),  USS  TRIPOLI  (LPH- 
10)  and  USS  PRINCETON  (CG-59)  all  struck  mines  and  received  severe  damage.  The 
costs  to  repair  the  SAMUEL  B.  ROBERTS  and  the  PRINCETON  totaled  over  $125 
million,  not  to  mention  the  loss  of  human  life.  This  staggering  figure  compared  to  the 
cost  of  the  mines  that  caused  the  damage  highlights  the  cost  effectiveness  of  mines.  Most 
recently,  during  Operation  Iraqi  Freedom  in  March  of  2003,  three  tugboats  were  stopped 
just  south  of  the  Iraqi  port  city  of  Umm  Qasr.  The  vessels  were  acting  suspiciously,  and 
were  found  to  be  carrying  70  mines  of  different  types.  The  tug  boats  and  the  barge  they 
were  towing  were  altered  so  that  mines  could  be  deployed  through  the  hull  to  avoid 
detection.  Some  of  the  mines  were  moored  mines  that  detonate  on  contact,  and  some 
were  bottom  influence  mines  that  are  detonated  by  a  ship’s  magnetic  or  acoustic 
signature. 
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Figure  1.  U.S.  troops  inspecting  mines  found  aboard  Iraqi  tug  boats. 


In  the  reverberation  limiting  acoustic  environment  of  the  littoral  regions  of  the 
world  (water  depths  less  than  100  meters),  mine  detection  can  prove  to  be  quite  difficult. 
As  proximity  of  the  mine  to  the  bottom  increases,  probability  of  detection  decreases  due 
to  complex  bathymetry  and  diverse  sediment  types.  Further  complicating  the  issue  are 
“bottom  boundary  layers”,  sometimes  called  “nepheloid  layers”  or  “suspended  sediment 
layers”.  These  layers  consist  of  suspended  sediments  that  occupy  the  lower  water  column, 
sometimes  up  to  10  m.,  and  lasting  for  several  weeks.  Their  effects  on  acoustic  detection 
are  undeniable.  Bottom  boundary  layers  can  make  a  mine  undetectable  and  increase  the 
threat  to  personnel  and  equipment.  Understanding  exactly  what  the  acoustic  effects  of 
bottom  boundary  layers  are  can  lead  to  the  development  of  acoustic  sensors  that  can  see 
through  to  the  sea  floor  and  aid  in  the  detection  of  the  potential  hazards  that  lie  within. 

In  this  thesis,  the  effects  of  a  suspended  bottom  sediment  layer  on  high  frequency 
acoustic  imagery  is  studied  using  the  Navy’s  Comprehensive  Acoustic  Simulation 
System  with  the  Gaussian  Ray  Bundle  (CASS/GRAB).  A  side  scan  sonar  image  of  a 
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mine-like  object  on  a  silty  clay  bottom  has  been  obtained  along  with  the  sound  velocity 
profile  and  all  bathymetric  data  corresponding  to  the  image. 

The  CASS/GRAB  model  is  run  with  the  appropriate  input  file  that  incorporates 
the  sound  velocity  profile,  bottom  type  and  bathymetric  data  and  the  scattering 
characteristics  of  a  mine-like  object.  The  bottom  synthetic  reverberation  field  calculated 
using  CASS/GRAB  is  taken  as  the  bottom  reverberation  as  well  as  the  echoes  from  the 
mine-like  object.  The  volume  scattering  strength  is  used  to  model  the  reverberation 
caused  by  the  suspended  sediment  boundary  layer.  The  CASS/GRAB  is  run  repeatedly 
while  changing  the  value  of  the  volume  scattering  strength  until  the  modeled  mine-like 
object  disappears,  which  leads  to  a  critical  value  representing  the  occurrence  of  the 
bottom  suspended  sediment  layer,  rendering  the  mine-like  object  undetectable.  The 
critical  values  of  volume  scattering  strength  are  then  taken  as  a  metric  of  the  acoustic 
impact  of  the  suspended  bottom  sediment  layer. 

The  remainder  of  the  thesis  is  organized  as  follows:  Chapter  II  discusses  the 
effects  that  suspended  sediment  has  on  acoustic  propagation.  Chapters  III  and  IV  describe 
the  CASS/GRAB  model  and  the  input  files  that  are  required  to  run  the  model  effectively. 
Chapters  V,  VI,  and  VII  discuss  the  results  of  this  work. 
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II.  WORKING  HYPOTHESIS  AND  PROCEDURE 


A,  SUSPENDED  SEDIMENT  AND  VOLUME  REVERBERATION 

One  version  of  the  active  sonar  equation,  where  there  is  zero  signal  excess,  is 
given  by 

SL-2TL  +  TS  =  RL  +  DT  (1) 

where  SL  represents  the  source  level  of  the  sonar,  TL  represents  transmission  loss,  TS 
represents  target  strength,  RL  represents  reverberation  level  and  DT  is  the  detection 
threshold  of  the  sonar.  Since  no  ambient  noise  information  was  available,  the  ambient 
noise  level  is  assumed  to  be  negligible.  Transmission  loss  is  multiplied  by  2  because  the 
radiated  sound  has  to  strike  the  target  and  then  return,  essentially  making  a  round  trip.  In 
the  detection  of  an  object  in  the  water,  the  source  level  and  target  strength  must  be  large 
enough  to  overcome  the  losses  of  energy  to  transmission  loss,  the  interference  due  to  the 
reverberation,  the  ability  of  the  sonar  processor  to  enhance  the  signal-to-noise  ratio  and 
the  ability  of  the  operator  to  recognize  objects. 

Transmission  loss  is  due  to  spreading  and  attenuation,  and  is  a  result  of  acoustic 
energy  traveling  through  water  (and  possibly  the  ocean  bottom)  and  reflecting  from  the 
ocean  bottom  and  surface.  There  are  three  types  of  reverberation;  volume  reverberation, 
sea-surface  reverberation  and  bottom  reverberation.  The  most  important  factor  that 
determines  reverberation  is  scattering  strength.  (Urick,  1983)  Scattering  is  the  reflection 
of  emitted  sound  from  material  in  the  water  column,  ocean  bottom  and  ocean  surface. 
The  mechanisms  that  cause  scattering  can  be  biological,  chemical  or  physical.  The 
governing  equation  for  volume  scattering  is: 

=101og^  (2) 

^ inc 

where  Sv  represents  the  scattering  coefficient,  in  decibels,  feat  represents  the  scattered 
intensity  of  the  sound  by  the  unit  volume,  and  Tnc  represents  the  incident  intensity  of  the 
sound.  When  there  are  particles  suspended  in  the  water  column,  the  scattering  also 
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depends  on  the  size  and  shape  of  the  seatterers.  In  this  work,  the  partieles  were  assumed 
to  be  spherieal  whieh  yields: 

s^  =  mog^  (3) 

4;r 

where  Sv  still  represents  the  seattering  eoeffieient  and  my  represents  the  baekseattering 
eross  seetion  of  a  unit  volume.  The  sum  of  all  seattered  energy  in  a  region  is  defined  as 
the  reverberation  (Uriek,  1983).  Not  only  does  the  size  and  shape  of  seatterers  affect  the 
reverberation,  but  the  concentration  of  seatterers  is  also  very  important.  The  more 
seatterers  there  are  in  a  unit  volume,  the  more  scattering  will  occur. 

Volume  reverberation  is  a  function  of  the  physical  properties  of  the  water  column 
including  any  sediment,  or  biological  material  that  may  be  suspended  there.  In  a 
suspended  bottom  sediment  layer,  large  quantities  of  sediment  remain  in  the  water 
column  for  up  to  several  weeks.  This  layer  significantly  affects  the  acoustic 
characteristics  of  the  water  that  it  is  suspended  in.  The  thicker  and  denser  a  suspended 
sediment  layer  becomes,  the  harder  it  would  be  for  a  sonar  to  penetrate  the  layer.  Another 
possible  impact  of  a  suspended  bottom  sediment  layer  would  be  changes  in  the 
temperature,  density  or  salinity  of  the  lower  water  column.  These  quantities  would  in  turn 
alter  the  sound  velocity  profile  and  might  prevent  acoustic  energy  from  reaching  a 
possible  mine.  The  focus  of  this  work  is  on  the  effects  of  a  suspended  sediment  layer  on 
the  volume  scattering  strength  component  of  reverberation. 


B,  PROCEDURE  TO  SIMULATE  SUSPENDED  SEDIMENT  LAYER 


The  working  hypothesis  of  this  thesis  is  “the  acoustic  impact  of  a  suspended 
sediment  layer  can  be  represented  by  high  values  of  volume  scattering  in  the  lower  water 
column”.  Using  an  actual  side-scan  sonar  image  as  a  reference,  an  acoustic  image  can  be 
created  using  CASS  and  the  bottom  can  be  effectively  “stirred-up”  to  create  a  suspended 
sediment  layer  that  did  not  previously  exist  in  the  image.  The  properties  of  this  simulated 
layer  can  then  be  changed  so  the  impact  of  the  layer  can  be  determined. 
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Figure  2  shows  a  representation  of  the  procedure  that  was  followed  to  recreate  a 
side-scan  sonar  image  using  CASS.  The  bathymetry  and  sound  velocity  profile  inputs  to 
CASS  were  taken  directly  from  the  side-scan  sonar  image.  Other  inputs,  including 
surface/bottom  parameters  and  volume  scattering  were  modeled  to  reflect  real  world 
conditions. 


Figure  2.  Steps  taken  to  create  a  total  reverberation  image  from  CASS,  compared 
to  the  side  scan  sonar  image. 

The  total  reverberation  image  is  compared  to  the  side-scan  sonar  image  for 
accuracy  before  the  suspended  bottom  sediment  layer  is  inserted.  Once  a  good 
comparison  has  been  made,  the  volume  scattering  in  the  lower  water  column  is  increased 
to  reflect  the  presence  of  the  suspended  sediment  layer. 

Figure  3  shows  the  decision  making  process  that  was  used  to  determine  if  the 
chosen  values  of  volume  scattering  were  adequate.  If  the  simulated  mine  object  is  still 
visible  after  increasing  the  volume  scattering,  then  the  simulated  sediment  layer  is  not 
strong  enough  to  represent  a  layer  that  would  hide  a  mine  from  a  sonar.  So  the  volume 
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scattering  was  increased  further.  This  procedure  was  performed  several  times  until  the 
mine  object  was  no  longer  visible. 


Figure  3.  Procedure  to  determine  whether  object  is  obscured  by  suspended 

sediment  layer 
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III.  COMPREHENSIVE  ACOUSTIC  SIMULATION  SYSTEM/  GAUSSIAN 

RAY  ACOUSTIC  BUNDLE 

A,  CASS  MODEL  FUNDAMENTALS 

CASS/GRAB  is  one  of  the  Navy’s  range-dependent  propagation,  reverberation, 
and  signal  excess  acoustic  models  that  can  be  used  in  both  the  active  and  passive  sonar 
equations.  The  CASS  model  acts  as  the  outer  shell  that  computes  range-independent 
monostatic  and  bistatic  active  signal  excess  calculations  and  calls  on  the  GRAB  subset 
for  the  computation  of  eigenrays  and  propagation  loss  (Figure  4). 


CASS 

Comprehensive  Acoustic 
System  Simulation 

Propagation  Model  1 :  FAME 


Propagation  Model  2:  GRAB 
Gaussian  Ray  Bundle 

Environmental  Interpolations 
Environmental  Model  Interpolations 
Surface  and  Bottom  Forward  Loss 
Volume  Attenuation 
Sound  Speed  Algorithms 


Propagation  Model  3:  COLOSSUS 
Propagation  Model  4:  AMOS  equations 
Backscatter  Models 
Reverberation 
Noise  Models 
Signal  to  Noise 
Signal  Excess 
Graphic  Displays 

System  Parameters  (Beamforming) 


Figure  4.  CASS/GRAB  Overview  (From  Keenan  et  al.2000) 
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In  the  GRAB  model,  travel  time,  source  angle,  target  angle,  and  phase  of  the  ray 
bundles  are  similar  to  those  values  for  a  classic  ray  path.  GRAB,  however,  uses  ray 
bundles  vice  a  ray  path.  This  allows  the  amplitude  of  the  Gaussian  ray  bundle  to  affect  all 
depths  to  some  degree,  whereas  with  classic  ray  paths,  amplitudes  are  local. 


B,  BATCH  FILE 

A  template  for  the  input  file  that  CASS  requires  to  compute  reverberation  can  be 
found  in  the  files  that  come  with  the  program.  A  more  specific  sample  of  the  input  file 
that  was  used  in  this  work  is  included  in  Appendix  A.  Once  created,  the  input  fide  must  be 
called  “input.dat”  for  CASS  to  read  it.  This  presented  a  slight  problem  since  the 
bathymetric  data  was  in  41  different  “lines”.  Each  line  represented  one  value  of  distance 
in  the  along  track  direction  and  had  several  values  of  cross  track  distance  and  depth 
associated  with  it.  Hence,  41  different  input  files  were  created  named  “inputl.dat” 
through  “input41.dat”. 

The  output  for  the  CASS  input  file  is  user  specified.  In  this  case,  the  output  file 
was  named  “revout.dat”  for  “reverberation  output”.  The  output  file  was  also  the  same  for 
each  file  and  needed  to  be  changed  to  something  that  made  data  retrieval  more  efficient. 

A  batch  file  was  created  that  took  “inputl.daf’  and  renamed  it  “input.dat”.  Then 
CASS  was  run  on  “input.dat”.  The  output  file  “revout.dat”  was  changed  to  “outputl.dat”. 
This  batch  file  continued  to  run  each  file  until  all  41  input  files  have  a  corresponding 
output  file.  This  batch  file  was  useful  because  it  allowed  for  multiple  runs  of  CASS  with 
varying  input  files  in  a  relatively  short  amount  of  time. 
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IV.  INPUT  PARAMETERS 


A,  ENVIRONMENTAL  PARAMETERS 

The  environmental  parameters  for  the  CASS  input  file  eonsisted  primarily  of  data 
taken  at  the  time  of  the  side-sean  sonar  image.  This  image  shows  an  objeet  that,  for  the 
purposes  of  this  work,  represents  a  “mine-like”  objeet.  The  objeet  is  assumed  to  be 
hollow  steel  and  sits  at  30  m  above  the  lower  boundary  of  the  image  and  27  meters  to  the 
left  of  the  right  boundary  of  the  image.  The  extent  of  the  image  is  approximately  60 
meters  in  the  y-direetion  and  50  meters  in  the  x-direetion.  The  following  seetions  will 
address  eaeh  part  of  the  environmental  parameters  one  at  a  time. 

1.  Bottom  Depth  Table 

The  values  for  the  bottom  depth  were  measured  at  the  time  the  image  was  taken. 
After  interpolating  between  adjaeent  values  to  improve  the  resolution,  values  were 
direetly  entered  into  the  table.  Values  of  depth  ranged  from  95  meters  to  77  meters.  A 
plot  of  bathymetry  eontours  is  ineluded  later  in  this  work. 

2,  Sound  Speed  Table 

The  values  for  the  sound  veloeity  profile  were  also  measured  at  the  time  of  the 
image.  Figure  5  shows  a  plot  of  the  sound  veloeity  vs.  depth. 
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Figure  5.  Sound  Velocity  Profile  for  region  above“mine-like”  object. 

3,  Bottom  Parameters 

Appendix  B  shows  that  the  grain  size  index  for  a  silty  clay  bottom  is  8.  APL/UW 
was  the  built-in  model  that  was  used  to  model  bottom  scattering.  Bottom  reflection 
effects  were  modeled  using  the  Rayleigh  scattering  model. 

4,  Surface  Parameters 

Wind  speed  at  the  time  the  image  was  created  was  not  recorded,  so  a  nominal 
wind  speed  of  1  m/s  was  reasonable  as  no  indication  of  significant  winds  were  evidenced 
in  the  side  scan  image.  APL/UW  was  the  built  in  model  that  was  chosen  to  model  surface 
reflection  and  the  surface  scattering. 

5,  Volume  Scattering  Strength 

For  this  work,  the  best  way  to  model  volume  scattering  strength  was  by  a  table 
that  was  directly  entered  into  the  input  file.  The  water  column  was  assumed  to  be 
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relatively  clear  above  77  meters,  with  an  initial  suspended  sediment  layer  characterized 
by  a  slightly  stronger  scattering  strength  below  78  meters. 


VOLUME  SCATTERING  STRENGTH  TABEE 


M 

DB/M 

0.00 

-95.00 

77.00 

-95.00 

78.00 

-65.00 

95.00 

-65.00 

Table  1.  Scattering  strength  table  for  CASS  input  file 

6.  Volume  Attenuation 

Out  of  all  the  models  available  in  CASS,  Erancois-Garrison  was  the  model  chosen 
for  volume  attenuation.  This  is  due  to  the  accuracy  that  this  model  handles  at  higher 
frequencies. 

B,  SONAR  PARAMETERS 

The  sonar  parameters  for  this  work  were  entered  to  reflect  that  of  a  high 
frequency  side  scan  sonar  at  100  kHz.  The  source  level  was  240  dB.  The  sonar  was  towed 
at  a  depth  of  30.4  meters  below  the  surface  and  had  a  pulse  length  of  0.001  seconds. 


C.  MODEL  PARAMETERS 

The  time  increment  for  modeling  should  not  exceed  one  half  of  the  pulse  length  to 
achieve  proper  resolution  of  each  time  step.  Since  the  total  distance  traveled  from  the 
sonar  to  the  end  of  the  image  is  approximately  50m,  the  total  reverberation  time  is  only 
0.12  seconds.  The  maximum  number  of  bottom  and  surface  reflections  modeled  was  set 
at  30  to  allow  for  some  interference  by  reflected  eigenrays. 
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V. 


CASS  REPRESENTATION  OF  SONAR  IMAGE 


The  first  task  that  was  to  be  aeeomplished  was  to  represent  the  side-sean  sonar 
image  (Figure  6)  through  reverberation  output  from  CASS. 


Figure  6.  Klein  5000  image  of  “mine-like”  objeet  on  silty  elay  bottom. 

Using  the  input  quantities  from  the  previous  seetion,  the  bateh  file,  and  a  Matlab© 
m-file  that  searehes  through  the  output  file  and  extraets  the  total  reverberation  data 
(Appendix  C.),  the  following  image  was  ereated  that  does  not  inelude  the  mine-like 
object; 
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Figure  7.  Reverberation  plot  of  bottom  without  mine  object  inserted.  X-axis  is 
cross  track  indices  represented  by  time  and  Y-axis  is  along  track  indices 
represented  by  distance. 

The  above  plot  shows  that  after  the  sound  leaves  the  sonar,  there  is  a  brief  time 
when  no  return  has  been  received.  The  first  blue  shading  is  the  return  from  volume 
reverberation.  The  first  red  line  and  the  yellow  field  that  follows  is  the  return  from  the 
surface  reverberation.  The  second  red  line  is  the  initial  return  from  the  bottom  (not  unlike 
a  fathometer).  The  final  feature  of  note  is  the  blue  “trench”  that  appears  near  the  far  edge 
of  the  plot.  This  is  because  the  sonar  was  located  at  a  depth  of  30.4  meters  which  was  just 
above  the  thermocline.  (Figure  5) 
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Now  that  there  is  an  accurate  depiction  of  the  bottom,  the  next  step  is  to  insert  the 
“mine-like”  object.  This  is  done  in  two  steps.  First,  the  bathymetry  must  be  changed  to 
reflect  a  mine  that  protrudes  above  the  bottom.  The  water  depth  in  the  vicinity  of  the 
object  is  87  meters.  The  size  of  the  object  in  the  image  is  5  meters  long,  3  meters  wide 
and  2  meters  high.  Therefore,  the  depth  in  the  vicinity  of  the  mine  will  change  to  85 
meters.  The  width  and  length  of  the  object  will  be  exaggerated  slightly  to  ensure  that  the 
object  is  hit  with  as  many  eigenrays  as  possible.  (Figures  8  and  9) 
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FIGURE  8.  Bathymetry  plot  without  mine  object. 
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FIGURE  9.  Bathymetry  plot  with  mine  object. 


The  second  step  to  entering  the  synthetic  “mine-like”  object  is  to  change  the 
bottom  parameters  to  reflect  hollow  steel  instead  of  silty  clay.  This  is  done  by  means  of 
inserting  a  bottom  scattering  strength  table.  (Table  1)  This  shows  for  all  grazing  angles, 
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the  scattering  strength  of  the  object  is  -1.13  dB/m  .  With  no  object  inserted  into  the 
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bathymetry,  there  was  only  one  bottom  parameter  for  the  input  file.  To  insert  an  object  in 
the  horizontal,  there  must  be  three  “environments”.  The  first  environment  is  the  same  as 
in  the  CASS  run  without  the  object.  Files  “outputl”  through  “outputlV”  cover  this 
regime.  The  second  “environment”  is  from  files  “outputlV”  thorough  “output22”.  This 
corresponds  to  an  object  that  is  roughly  8  meters  long  in  the  y-direction.  In  the  x- 
direction,  the  altered  bottom  scattering  values  occupy  values  from  27.8  meters  to  32.8 
meters,  producing  an  image  that  is  5  meters  wide.  The  original  object  was  5  meters  long 
by  3  meters  wide.  This  simulated  image  is  slightly  larger,  but  still  representative  of  the 
original. 

Now,  CASS  is  run  again  with  the  altered  files.  The  following  image  is  the  product 
of  this  attempt; 
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FIGURE  10.  Reverberation  plot  of  bottom  with  mine  object  inserted.  X-axis  is 
cross  track  indices  represented  by  time  and  Y-axis  is  along  track  indices  represented  by 

distance 
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Clearly,  the  object  is  visible  in  the  reverberation  imagery.  It  is  located  at  x- 
cooordinate  150,  y-coordinate  23.  All  of  features  that  were  in  the  image  without  the  mine 
are  still  present,  but  now  the  object  itself  and  even  an  echo  from  the  object  are  visible. 
The  echo  is  the  orange  shaded  area  located  near  the  left  edge  of  Figure  10.  The  value  of 
this  image  is  that  now  there  exists  a  synthetic  replica  of  a  real  image  that  inputs  that  can 
be  changed  to  examine  the  effects  of  various  parameters. 
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VI. 


ACOUSTIC  IMPACT  OF  SUSPENDED  SEDIMENT  LAYER 


As  shown  in  the  input  parameters  seetion,  the  volume  seattering  table  showed  a 
slightly  elevated  seattering  value  for  the  depths  below  78  meters  (Table  1).  This  elevation 
was  part  of  a  typieal  volume  seattering  profile.  Sediment  in  the  water  column  would 
increase  the  volume  scattering  and  ultimately  the  volume  reverberation.  Since  the 
simulated  side-scan  sonar  image  with  an  object  inserted  has  inputs  that  can  be  changed,  a 
critical  value  can  be  determined  as  to  how  much  volume  scattering  in  a  simulated 
suspended  bottom  sediment  layer  will  render  the  “mine-like”  object  undetectable. 
Volume  attenuation  and  changes  in  the  sound  velocity  profile  will  also  have  an  effect,  but 
they  are  not  addressed  in  this  work. 

Starting  with  the  input  files  that  had  the  object  in  place,  the  initial  value  of  -65 
dB/m  was  changed  to  -60  dB/m  .  After  plotting  the  output,  little  change  was  noted.  For 
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the  next  attempt,  the  volume  scattering  was  changed  to  -50  dB/m  .  Changes  to  the  plot 
began  to  appear.  The  reverberation  values  around  the  object,  in  the  vicinity  of  where  the 
bottom  reverberation  appeared  before,  were  increasing.  This  was  significant  progress,  but 
the  object  was  still  visible.  Several  more  attempts  were  made  until  at  -30  dB/m  ,  the 
object  was  almost  gone.  (Figure  1 1) 
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FIGURE  11.  Reverberation  plot  of  bottom  with  mine  object  inserted.  X-axis  is 
cross  track  indices  represented  by  time  and  Y-axis  is  along  track  indices  represented  by 

-5 

distance  with  mine  object  inserted  and  suspended  sediment  layer  with  -30  dB/m 
volume  scattering  strength. 

The  object  is  still  visible  in  Figure  11,  but  it  is  beginning  to  approach  that 
threshold  where  mine  detection  equipment  might  not  be  able  to  distinguish  the  object 
from  the  surrounding  bathymetry.  Still,  the  goal  of  the  experiment  was  to  find  the  values 
of  volume  scattering  that  rendered  the  object  undetectable.  Next,  volume  scattering 
values  were  changed  and  1  dB/m  intervals  until  finally  at  -22  dB/m  ,  the  object  was 
completely  obscured  by  the  suspended  sediment  layer.  (Figure  12) 
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Reverberation  (dB)  with  object  present  at  -22  dB/m 
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FIGURE  12.  Reverberation  plot  of  bottom  with  mine  object  inserted.  X-axis  is 
cross  track  indices  represented  by  time  and  Y-axis  is  along  track  indices  represented  by 
distance  with  mine  object  inserted  and  suspended  sediment  layer  with  -22  dB/m^  volume 
scattering  strength 
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VII.  CONCLUSIONS 


It  has  been  shown  that  the  CASS/GRAB  model  ean  be  used  to  ereate  a  side  sean 
sonar  image  by  mapping  the  modeled  reverberation  and  that  a  synthetie  mine-like  objeet 
ean  be  inserted  into  the  input  file  to  refieet  a  real  objeet.  Also,  for  this  partieular  instanee, 
volume  seattering  values  of  -30  dB/m  to  -22  dB/m  render  the  objeet  aeoustieally 
undeteetable.  These  values  are  speeifie  to  this  ease  only  and  do  not  represent  universal 
values.  The  results  of  this  study  must  also  be  tempered  with  the  knowledge  that  ehanges 
in  volume  attenuation  and  sound  speed,  eorresponding  to  the  ehange  in  volume 
seattering,  were  not  ineluded. 

This  researeh  was  valuable  beeause  it  provides  a  produet  that  ean  model  the 
effeets  of  a  suspended  sediment  layer  on  side  sean  sonar  imagery.  Given  the  appropriate 
inputs,  this  produet  ean  provide  results  for  a  taotieally  signifieant  issue  for  the  mine 
warfare  eommunity.  While  the  development  of  this  produet  is  signifieant,  several 
shortfalls  remain.  First,  the  proeess  by  whieh  the  environment  and  the  objeet  are  modeled 
is  eumbersome.  Seeond,  the  appropriate  volume  attenuation  and  sound  speed  must  also 
be  used  with  into  this  produet.  Follow-on  efforts  should  provide  solutions  to  the  two 
issues  listed  above.  Equally  (or  more  important)  is  a  thorough  study  of  the  relationship 
between  the  suspended  sediment  layer  density  and  type  (e.g.,  sand,  silt  or  clay),  particle 
density  in  the  layer,  associated  volume  scattering  and  attenuation,  and  changes  in  the 
sound  speed  profile. 
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APPENDIX  A.  CASS/GRAB  MODEL  INPUT  CARD  FOR  IMAGE  WITHOUT 
MINE-LIKE  OBJECT  OR  SUSPENDED  SEDIMENT  LAYER 


Direct  Output 


OUTPUT  FILE  =  REVOUT 
RESET  OUTPUT  DEVICE 

VERTICAL  ANGLE  MINIMUM  =  0  DEG 
VERTICAE  ANGEE  MAXIMUM  =  90  DEG 
VERTICAL  ANGLE  INCREMENT  =  I  DEG 


Erequency 


EREQUENCY  MINIMUM  =  100000.00  HZ 
EREQUENCY  MAXIMUM  =  100000.00  HZ 
BANDWIDTH  TABEE  =  1 .00  HZ 


BATHYMETRY 


BOTTOM  DEPTH  TABEE 
M  M  M 

0  0  84.94 

1.284  0  85.5 

2.568  0  86.06 

3.852  0  86.34 
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5.136 

0 

86.62 

6.345 

0 

86.905 

7.554 

0 

87.19 

8.838 

0 

87.095 

10.122 

0 

87 

11.406 

0 

87.28 

12.691 

0 

87.56 

13.899 

0 

87.845 

15.105 

0 

88.13 

16.392 

0 

88.41 

17.676 

0 

88.69 

18.960 

0 

89.065 

20.245 

0 

89.44 

21.453 

0 

89.72 

22.662 

0 

90 

23.946 

0 

90.47 

25.230 

0 

90.94 

26.515 

0 

91.41 

27.799 

0 

91.88 

29.008 

0 

91.97 

30.216 

0 

92.06 

31.500 

0 

91.5 

32.785 

0 

90.94 

34.069 

0 

91.41 

35.353 

0 

91.88 

36.562 

0 

91.505 

37.770 

0 

91.13 

39.055 

0 

92.16 

40.339 

0 

93.19 

41.623 

0 

93.375 

42.907 

0 

93.56 

44.116 

0 

94.31 
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45.325 

0 

95.06 

46.609 

0 

94.97 

47.893 

0 

94.88 

49.177 

0 

95.065 

50.461 

EOT 

0 

95.25 

Sound  speed  Table 


ENVIRONMENT  X  COORDINATE  =  0.0  M 
ENVIRONMENT  Y  COORDINATE  =  0.0  M 


SOUND  SPEED  TABEE 

M 

M/S 

0.00 

I5I9.998I 

5.00 

I5I9.998I 

6.00 

1520.0421 

7.00 

1520.0862 

8.00 

1520.1026 

9.00 

I520.II90 

10.00 

1520.1076 

1 1. 00 

I520.I5I7 

12.00 

1520.1958 

13.00 

1520.2233 

14.00 

1520.2285 

15.00 

1520.2172 

16.00 

1520.2336 

17.00 

1520.2500 

18.00 

1520.2664 

19.00 

1520.2827 

20.00 

1520.3103 

29 


21.00 

1520.3543 

22.00 

1520.3707 

23.00 

1520.3871 

24.00 

1520.4035 

25.00 

1520.4475 

26.00 

1520.4639 

27.00 

1520.4803 

28.00 

1520.5078 

29.00 

1520.5242 

30.00 

1520.5683 

31.00 

1520.5570 

32.00 

1520.5734 

33.00 

1520.6174 

34.00 

1520.6338 

35.00 

1520.7830 

36.00 

1521.0372 

37.00 

1521.0483 

38.00 

1520.8550 

39.00 

1520.4563 

40.00 

1519.8601 

41.00 

1519.2093 

42.00 

1516.5164 

43.00 

1515.6759 

44.00 

1515.3182 

45.00 

1515.3036 

46.00 

1515.2825 

47.00 

1514.8935 

48.00 

1512.9059 

49.00 

1512.3495 

50.00 

1512.0346 

51.00 

1511.8022 

52.00 

1511.1448 

30 


53.00 

1511.0015 

54.00 

1510.8884 

55.00 

1510.6417 

56.00 

1510.5972 

57.00 

1510.5563 

58.00 

1510.6303 

59.00 

1510.6662 

60.00 

1510.6980 

61.00 

1510.6953 

62.00 

1510.7117 

63.00 

1510.7281 

64.00 

1510.7139 

65.00 

1510.6805 

66.00 

1510.6662 

67.00 

1510.6327 

68.00 

1510.5377 

69.00 

1510.5234 

70.00 

1510.5206 

71.00 

1510.5063 

72.00 

1510.4919 

73.00 

1510.4467 

74.00 

1510.4324 

75.00 

1510.4180 

76.00 

1510.4151 

77.00 

1510.4316 

78.00 

1510.4480 

79.00 

1510.4027 

80.00 

1510.3883 

81.00 

1510.3739 

82.00 

1510.3903 

83.00 

1510.4183 

84.00 

1510.4039 
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85.00 

1510.3701 

86.00 

1510.3865 

87.00 

1510.4029 

88.00 

1510.3885 

89.00 

1510.3856 

90.00 

1510.3402 

91.00 

1510.2561 

92.00 

1510.2416 

93.00 

1510.2580 

94.00 

1510.2744 

95.00 

1510.2909 

96.00 

1510.2763 

97.00 

1510.2928 

98.00 

1510.2163 

99.00 

1510.2017 

100.00 

1510.2182 

101.00 

1510.2036 

102.00 

1510.1890 

103.00 

1510.2054 

104.00 

1510.1599 

105.00 

1509.9590 

106.00 

1509.7888 

107.00 

1509.6496 

108.00 

1509.6348 

109.00 

1509.5577 

110.00 

1509.3868 

111.00 

1509.1334 

112.00 

1508.8560 

113.00 

1508.4836 

114.00 

1508.3856 

115.00 

1508.2643 

116.00 

1507.8065 
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117.00 

1507.4740 

118.00 

1507.3517 

119.00 

I507.I25I 

120.00 

1506.4675 

121.00 

1505.7128 

122.00 

1504.5219 

123.00 

1498.3821 

124.00 

1490.8919 

125.00 

1485.0349 

126.00 

1482.3583 

127.00 

1480.0970 

128.00 

1478.3381 

EOT 


Surface  Reflection  Coefficient 


WIND  SPEED  =  2.00  KNOTS 

SUREACE  REEEECTION  COEEEICIENT  MODEE  =  APE/UW 

EUNCTION  SYMBOE  =  SRE  RFE 

FUNCTION  UNIT  =  DB 

PRINT  FUNCTION  VS  VERTICAE  ANGEE 


Volume  Attenuation 


VOEUME  ATTENUATION  MODEE  =  FRANCOIS 
FUNCTION  SYMBOE  =  VEM  ATN 
FUNCTION  UNIT  =  DB 
PRINT  FUNCTION  VS  VERTICAE  ANGEE 
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Surface  Scattering 


SURFACE  SCATTERING  STRENGTH  MODEE  =  APE/UW 

EUNCTION  SYMBOE  =  SRE  STR 

EUNCTION  UNIT  =  DB 

PRINT  FUNCTION  VS  VERTICAE  ANGEE 


Volume  Scattering 


VOEUME  SCATTERING  STRENGTH  TABEE 
M  DB//M 
0.00  -95.00 
77.00  -95.00 
78.00  -65.00 
95.00  -65.00 
EOT 


Bottom  Parameters 


BOTTOM  SEDIMENT  GRAIN  SIZE  INDEX  =  8 

BOTTOM  SCATTERING  STRENGTH  MODEE  =  APE/UW 

EUNCTION  SYMBOE  =  BTM  REE 

EUNCTION  UNIT  =  DB 

PRINT  FUNCTION  VS  VERTICAL  ANGLE 

BOTTOM  REELECTION  COEEEICIENT  MODEL  =  RAYLEIGH 
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FUNCTION  SYMBOL  =  BTM  STR 

FUNCTION  UNIT  =  DB 

PRINT  FUNCTION  VS  VERTICAL  ANGLE 


Source  parameters 


SOURCE  EEVEE  TABEE  =  240.00  DB 
PUESEEENGTH=  0.001  S 
TRANSMITTER  DEPTH  =  30.4  M 
SOURCE  DEPTH  =  30.4  M 


Transmitter  Beam  Pattern 


TRUE  TRANSMITTER  BEARING  =  0  DEG 
TRUE  TRANSMITTER  HEADING  =  0  DEG 
TRANSMITTER  HORIZONTAE  BEAMWIDTH  MODEE  =  TABEE 
TRANSMITTER  HORIZONTAL  BEAMWIDTH  TABEE 
KHZ  DEG 
100.00  .15 

EOT 

TRANSMITTER  BEAM  PATTERN  MODEE  =  TABEE 
TRANSMITTER  BEAM  PATTERN  TABEE 
DEG  DB 
-90.0  -900.00 
-86.0  -900.00 
-85.0  -20.00 
-28.0  -20.00 
-12.0  0.00 
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4.0  0.00 

20.0  0.00 
36.0  -20.00 
90.0  -20.00 
EOT 


Receiver  Beam  Pattern 


BEARING  ANGEE  MINIMUM  =  0  DEG 
BEARING  ANGEE  MAXIMUM  =  0  DEG 
BEARING  ANGEE  INCREMENT  =  I.O  DEG 
RECEIVER  HORIZONTAL  BEAMWIDTH  MODEL  =  TABLE 
RECEIVER  HORIZONTAL  BEAMWIDTH  TABLE 
KHZ  DEG 
100.00  .15 

EOT 

RECEIVER  BEAM  PATTERN  MODEL  =  TABLE 
RECEIVER  BEAM  PATTERN  TABLE 
DEG  DB 


-90.0 

-900.00 

-82.0 

-900.00 

-81.0 

-20.00 

0.0 

-20.00 

1.0 

0.00 

15.0 

0.00 

16.0 

-20.00 

90.0 

-20.00 

EOT 
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Model  and  Run  Parameters 


EIGENRAY  MODEE  =  GRAB 
MAXIMUM  BOTTOM  REEEECTIONS  =  30 
MAXIMUM  SURE  ACE  REEEECTIONS  =  30 
TIME  MINIMUM  =  0.0  S 

TIME  MAXIMUM  =  .  1 2  S 

TIME  INCREMENT  =  0.0005  S 

RANGE  MINIMUM  =  0.0  M 

RANGE  MAXIMUM  =  50.0  M 

RANGE  INCREMENT  =  0.375  M 
VERTICAE  ANGEE  MINIMUM  =  -89  DEG 
VERTICAE  ANGEE  MAXIMUM  =  89  DEG 
VERTICAE  ANGEE  INCREMENT  =  .  I  DEG 
ACTIVE  MODE  =  MONOSTATIC 


bottom  depth  eigenray  seetion 


TARGET  DEPTH  =  BOTTOM 
EIGENRAY  EIEE  =  EIGBTM 
EIGENRAY  ADDITION  =  RANDOM 
COMPUTE  EIGENRAYS 


surfaee  depth  eigenray  seetion 


TARGET  DEPTH  =  SURE  ACE 
EIGENRAY  EIEE  =  EIGSRE 
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EIGENRAY  ADDITION  =  RANDOM 
COMPUTE  EIGENRAYS 


Target  depth  eigenray  seetion 


TARGET  DEPTH  =  20.0  M 
SCATTERING  EAYER  THICKNESS  =  40.0  M 
EIGENRAY  EIEE  =  EIGOOI 
EIGENRAY  ADDITION  =  RANDOM 
COMPUTE  EIGENRAYS 
TARGET  DEPTH  =  60.0  M 
SCATTERING  EAYER  THICKNESS  =  40.0  M 
EIGENRAY  EIEE  =  EIG002 
EIGENRAY  ADDITION  =  RANDOM 
COMPUTE  EIGENRAYS 

TARGET  DEPTH  =  81  M 
SCATTERING  EAYER  THICKNESS  =  15.0  M 
EIGENRAY  EIEE  =  EIG003 
EIGENRAY  ADDITION  =  RANDOM 
COMPUTE  EIGENRAYS 


Compute  Reverberation 


EATHOMETER  RETURN  MODEE  =  NB 
REVERBERATION  MODEE  =  NB 
RESET  REVERBERATION 
COMPUTE  EATHOMETER  RETURNS 
EIGENRAY  EIEE  =  EIGSRE 
COMPUTE  SURE  ACE  REVERBERATION 
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EIGENRAY  EIEE  =  EIGBTM 
COMPUTE  BOTTOM  REVERBERATION 
EIGENRAY  EIEE  =  EIGOOl 
COMPUTE  VOEUME  REVERBERATION 
EIGENRAY  EIEE  =  EIG002 
COMPUTE  VOEUME  REVERBERATION 
EIGENRAY  EIEE  =  EIG003 
COMPUTE  VOEUME  REVERBERATION 

PRINT  REVERBERATION  VS  TIME 

END 
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APPENDIX  B,  BOTTOM  TYPE  GEO-ACOUSTIC  PROPERTIES  FROM 

NAVO  (1999) 


Bottom 

CompositKXt 

Gfain 

Size 

Long  Name 

Sound 

Speed 

Ratw 

Wave 

Number 

Ratio 

Volume 

Parameter 

Spectral 

Exponent 

Spectral 

Parameter 

BOaOER 

ROUGH  ROCK 

2  50C' 

2  5000' 

001374 

0  002C 

325 

0  20693C' 

RXK 

-70 

RXK 

2500 

2500C 

0  01374 

0  002C' 

325 

001862C' 

COBBLE 

2  500 

1  8000' 

001374 

0  0020 

325 

0  016000 

GRAVa 

Km 

GRAVa 

2500 

18000 

001374 

0  0020 

325 

0015000 

-30 

PEBBLE 

2  500 

1800C' 

001374 

0  002C 

325 

001600C' 

Km 

SANDY  GRAVEL 

2  492 

1  3370 

001705 

0  0020 

325 

0012937 

2  401 

1  3057 

0002Q 

325 

0  010573 

KOI 

MUDDY  SANDY  GRAVEL 

2314 

1277S 

001630- 

325 

0  0086G2 

mm 

COARSE  SAfO 

2  231 

12503 

001638 

00020 

325 

00*95“ 

05 

GRAVELLY  SAND 

2  231 

12503 

0  002C 

325 

000^ 

■Di 

GRAVELY  MUDDY  SAt© 

2151 

12241 

0.01645 

0  0020 

325 

0  00558’ 

SAND 

15 

SAt© 

1845 

11782 

001624 

0  0020 

325 

0  004446 

15 

KCDIUMSAND 

1845 

11782 

001624 

0  0020 

325 

0004446 

KM 

MUDDY  GRAVEL 

1615 

11396 

001610 

0  0020 

325 

0003496 

25 

FINE  SAND 

1  451 

11073 

001602 

00020 

325 

0002715 

SILTY  SAND 

1  451 

1 1073 

001602 

00020 

325 

0002715 

KM 

MUDDY  SAND 

1.339 

001728 

0  0020- 

325 

0002070 

35 

VERY  FINE  SAND 

1268 

1,0568 

001875 

00020 

325 

0001544 

■M 

CLAYEY  SAND 

1224 

10364 

002019 

00020 

325 

0001119 

45 

COARSE  SILT 

1195 

10179 

0  02158 

00020 

325 

0000781 

50 

SANDY  SILT 

1 169 

09999 

001261 

0  0020 

3.25 

0000518 

5,5 

H€DIUMSLT 

1 149 

09885 

00010 

325 

0000518 

SAND-SILT©LAY 

1  149 

09885 

000676 

00010 

325 

0000518 

SILT 

■M 

SILT 

1149 

000386 

00010 

325 

0000518 

SANDY  MUD 

1149 

09873 

0  00386 

00010 

325 

0  000518 

65 

Fl^e  SILT 

1.148 

09861 

0  00306 

0  0010 

325 

0000518 

65 

CLAYEY  SILT 

1 148 

09661 

0  00306 

00010 

325 

0  000518 

MUD 

mm 

SANDY  CLAY 

1 147 

09849 

000242 

00010 

325 

0.000518 

75 

VERY  FINE  SILT 

1 147 

09837 

000194 

00010 

3,25 

0,000518 

■M 

SILTY  CLAY 

1146 

09824 

000163 

325 

0000518 

CLAY 

KM 

CLAY 

1145 

09600 

0.00148 

00010 

325 

0000516 

■EEI 

1145 

KEiiM 

000148 

00010 

325 
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APPENDIX  C.  MATLAB©  M-FILE  TO  RETRIEVE  REVERBERATION 

INFORMATION 


%  file  =  PlotCASSReverb  all.m 
% 

clear 

defiow  =  0.0;  %  min  value  is  Zero 

for  kk  =  1:41 

fname  =  ['output'  int2str(kk)  '.DAT'] 
fid  =  fopen(  fname, 'r'); 

% 

%  Acquire  reverberation  data. 

% 

string  =  fgetl(fid); 
boogie  =  1; 
while  boogie  ==  1 
string  =  fgetl(fid); 
if  length(string)  >  25 
if  strmg(24:25)  ==  'NB' 
boogie  =  0; 
end 
end 
end 

string  =  fgetl(fid); 
string  =  fgetl(fid); 
string  =  fgetl(fid); 
string  =  fgetl(fid); 
for]  =  1:241 
string  =  fgetl(fid); 

A(l,j)  =  str2num(strmg(l:12));  %  Col  1  —  Time 
A(2,j)  =  max(deflow,str2num(strmg(61 :72)));  %  Col  6 


end 
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fclose(fid); 

Atime  =  A(l,:);  %  Min  time  is  1st  value  in  row  (Left) 

index  =  find(A(2,:)  ==  0);  %  Find  zeros 

A(2, index)  =  NaN;  %  replaee  zero  with  NaN 

%Atotal(kk,:)  =  A(2,:);  %Left-to-right  row  veetor  with 

A1  in  Row  1  (top  of  matrix) 

Atotal(42-kk,:)  =  A(2,:);  %Left-to-right  row  veetor 
with  A1  in  Row  41  (bottom  of 
matrix) 

End 
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